From nonlinear scaling fields to critical amplitudes 



Y. Meurice and S. Niermann 
Department of Physics and Astronomy, The University of Iowa, Iowa City, Iowa 52242, 

USA 



Abstract 



We propose to combine the nonlinear scaling fields associated with the 
high-temperature (HT) fixed point, with those associated with the unstable 
fixed point, in order to calculate the susceptibility and other thermodynamic 
quantities. The general strategy relies on simple linear relations between the 
HT scaling fields and the thermodynamic quantities, and the estimation of 
RG invariants formed out of the two sets of scaling fields. This estimation 
requires convergent expansions in overlapping domains. If, in addition, the 
initial values of the scaling fields associated with the unstable fixed point 
can be calculated from the temperature and the parameters appearing in 
the microscopic Hamiltonian, one can estimate the critical amplitudes. This 
strategy has been developed using Dyson's hierarchical model where all the 
steps can be approximately implemented with good accuracy. We show nu- 
merically that for this model (and a simplified version of it), the required 
overlap apparently occurs, allowing an accurate determination of the critical 
amplitudes. 

Keywords: Renormalization group, scaling fields, high-temperature expan- 
sion, hierarchical model, normal forms, critical amplitudes, crossover. 



I. INTRODUCTION, MOTIVATIONS AND MAIN RESULTS 



It is well-known that the magnetic susceptibility of a spin model near its critical tem- 
perature can be parametrized as 

X = (/5c - ^)-\Ao + -^)^ + ...). (1.1) 

In this expression, the exponents 7 and A are universal and can be obtained from the 
calculation of the eigenvalues of the linearized renormalization group (RG) transformation. 
On the other hand, the critical amplitudes Aq, Ai, . . . are functions of the microscopic 
details of the theory. One can find universal relations among these amplitudes and the 
ones associated with other thermodynamic quantities, however these relations do not fix 
completely the amplitudes. In the end, if we want a quantitative estimate of a particular 
amplitude, we need to perform a calculation which requires a knowledge of many details of 
the RG flows. Such a calculation is in general a difficult, nonlinear, multivariable problem. In 
this article we propose a general strategy based on the construction of nonlinear scaling fields 
associated with several fixed points, to calculate the critical amplitudes, and we demonstrate 
its feasibility in the case of Dyson's hierarchical model. 

A common strategy in problems involving nonlinear fiows near a singular point, is to 
construct a new system of coordinates for which the governing equations become linear. It 
seems intuitively clear that if the original problem is sufficiently nontrivial, normal form 
methods can only work in some limited way, locally, because the fiows of the nonlinear 
problem have global properties which do not match those of the linear fiows. A well-known 
argument for the inadequacy of such procedure (which extends beyond the special case of an 
expansion near a singular point) was provided by Poincare in the context of perturbed 
integrable Hamiltonians. He discovered that even though it is possible to write a formal per- 
turbative series for the action-angle variables, some coefficients have "small denominators" , 
and generically, the series are ill-defined. However, under some restrictions (formulated ac- 
cording to some appropriate version of the K. A. M. theorem 0), perturbation theory can 
still provide interesting information. 

Almost thirty years ago, Wegner [Q, introduced quantities that transformed multiplica- 
tively under a RG transformation. He called them "scaling fields" and we will use his 
terminology in the following. Sometimes, one also uses the terminology "nonlinear scaling 
field" to distinguish them from the linear ones (see section || for details). In the following, 
"scaling fields" mean the nonlinear ones and we will use the terminology "linear scaling 
fields" when necessary. These fields play a central role in the presentation of the basic ideas 
of the RG. They appear in almost any review on the subject (see for instance Ref. ||^). As 
in the case of Hamiltonian dynamics, there exists a formal series expansion for the scaling 
variables (see Eq. (4.9) in Ref. |^). It is commonly assumed that the functions defined with 
this procedure are analytic, at least within a certain neighborhood of the fixed point. How- 
ever, for most non-trivial models, it is very difficult to prove this assumption. In particular, 
it is difficult to address the question of small denominators because it requires an accurate 
calculation of the eigenvalues of the linearized RG transformation. 

If the small denominator problem can be controlled and if some local expansion is well- 
defined, there remain several important global issues. What is the domain of convergence of 
this expansion? How does the accuracy of an expansion with a finite number of terms evolve 
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when we move away from the fixed point? Can different expansions have overlapping do- 
main of convergence? These important global issues are rarely discussed because of practical 
limitations: in crossover regions, we need large order expansions in many variables. Unfor- 
tunately, this problem has to be faced if we want to calculate all the critical amplitudes. In 
this article, we propose a general strategy to calculate directly the critical amplitudes. This 
strategy has been developed using Dyson's hierarchical model, where large order expansions 
in many variables are practically feasible. All the numerical calculations presented hereafter 
were done with this model (or a simplified version of it). 

The general point of view that we want to advocate here is that one should combine 
different sets of scaling fields. Even though the scaling fields are almost always constructed in 
the vicinity of Wilson's fixed point, they can in principle be constructed near any other fixed 
point. If one can find some overlap among the domains of convergence of these expansions it 
is possible to reconstruct the fiows, given their initial values. In other words, we would like 
to develop a new analytical approach to complement the existing methods used to deal with 
the crossover between fixed points, namely, the Monte Carlo method a combination 

of field-theoretical methods and mean field calculations P,p!0| or the study of the entropy 



associated with the RG fiows 111 



In the following, we concentrate on the study of the RG fiows in the symmetric phase 
of spin models having a nontrivial unstable fixed point. Our goal is to calculate the critical 
amplitudes by constructing the scaling fields near the three relevant fixed points: the Gaus- 
sian fixed point (if relevant), the unstable fixed point (sometimes called the IR fixed point or 
Wilson's fixed point), and the high-temperature (HT) fixed point. The idea is represented 
schematically in Fig. |1|. 

We propose to follow three specific steps to achieve this goal. These steps correspond to 
a construction in backward order, starting with the fiows near the HT fixed point and ending 
with the initial conditions. First, we express the thermodynamic quantities in terms of the 
scaling fields of the HT fixed point. Second, we use the scaling fields of the unstable IR fixed 
point, to write the thermodynamic quantities as a main singularity times a RG invariant 
quantity constructed out the two scaling fields. Third, we calculate the initial values of 
the scaling fields associated with the unstable fixed point in terms of the basic parameters 
appearing the microscopic Hamiltonian. These three steps are explained in more detail in 
Section 0. This section is essential to understand the general ideas and the notations used 
later. It should be noted that the first two steps are independent of the initial conditions 
and are in some sense universal. On the other hand, the third step provides the initial data 
as a function of the basic parameters such as the temperature or the "bare parameters" 
of Landau- Ginzburg models. Consequently, the method used to implement the third step 
depends on the initial measure considered. For instance, if we start with an initial measure 
near the Gaussian fixed point, perturbative field theoretical methods (Feynman diagrams) 
will be used. 

Each of these three steps can in principle be implemented for spin models with nearest 
neighbor interactions in three dimensions, by using available expansions to describe the 
flows in each region. However, in order to discuss the behavior of the expansions in the 
crossover region, we need large order expansions. Despite the existence of increasingly 
sophisticated methods used for various expansions in nearest neighbor models (see e.g., 
Ref. 10), it is still a major time investment to perform these expansions. In order to 
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test the feasibility of the procedure, we have considered approximations for which large 
order expansions (see e.g., Refs. |l^,14l) can be reached more easily. Namely, we have used 



hierarchical approximations where only the local part of the measure is renormalized under 

or 



a RG transformation. Well-known examples are the "approximate recursion formula" [|15 
Dyson's hierarchical model (HM) |T^|T^. Despite this approximation, the nonlinear aspects 
of the problem are still nontrivial. To fix the ideas, the calculations that we have performed 
required the determination of several thousands of coefficients in various expansions. 



The relevant facts about the HM and its RG transformation are reviewed in Section [II 



Before embarking in a multivariable calculation, we have first considered a simplified version 
involving only one variable |18[ where the small denominator problem is obviously absent. 
This model discussed in Section is simply a quadratic map with two fixed points, one 
stable and one unstable. Series expansions for the scaling fields of this simplified model, and 
their inverse, can easily be constructed. Treating these numerical series with well-known 
estimators |]TU[, we obtain radii of convergence and exponents in very good agreement with 
what we can infer by using general arguments. These numerical results will be used as 
references when a similar analysis is conducted later for the HM. The most important result 
of Section |^ is illustrated Fig. |^ which shows that the expansions of the scaling fields scale 
accurately in overlapping regions. 

The rest of the article is devoted to the HM with one spin component per site (as in the 
Ising model). For this model, the existence of a non-trivial unstable fixed point has been 
proven rigorously [p0| , pl| . Significant results have been obtained regarding the local existence 



of scaling fields for the RG map near the unstable fixed point by Collet and Eckmann ||2T[ 



and Koch and Wittwer |]22|. In addition, Ref. contains a mathematical justification of 
the polynomial approximations that we have used to perform large order HT expansions 
|T^,|T^ or direct numerical calculations p5| , P^ . In the following, the parameter playing the 



role of the dimensionality (see Section |lTTD will be tuned in such way that a Gaussian massless 
field scales exactly as in three dimensions. In other words, we will work at an intermediate 
value between the upper and lower critical dimensions and the e-expansion is not obviously 
useful. For this particular choice of the dimensionality parameter, the numerical value of 
the unstable fixed point is known with great precision in a specific system of coordinates 



| 25| and the question of small denominators studied in Ref. p6| . 

Using these results we first provide an explicit construction of the two sets of scaling 
fields and show that the first two steps can be implemented. We found simple linear rela- 
tions between the derivatives of the free energy and the HT scaling fields. The RG invariants 
discussed above can be calculated accurately (see Fig. [M]) because the scaling fields associ- 
ated with the two fixed points scale accurately in overlapping regions as shown in Fig. |T^. 
The third step was performed in the case of an initial Ising measure where the temperature 
is the only free parameter. Putting everything together, we can calculate the leading and 
subleading amplitude of the susceptibihty and have very good agreement with previous nu- 
merical calculations p4|. This demonstrates the feasibility of the proposed startegy. In the 
conclusions, we discuss our plans to extend this method for perturbative initial conditions 
and for nearest neighbor models. 

A few words of caution. We want to make clear that the main result presented in this 
article is a calculation of the leading and subleading critical amplitudes of the magnetic 
susceptibility of the hierarchical model. This calculation shows that the strategy discussed 
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in the coming section II can be implemented successfully for this particular model. The 
fact that part of the introduction and Section II are written for a general spin model should 
not be interpreted as the statement that it is straightforward to implement the strategy for 
other models than the hierarchical model. The reader whose main interest is the hierarchical 
model may read section II as a general description of what is done in the rest of the paper. 
On the other hand, the reader interested in using the general strategy for other models, 
should be aware that we do not provide a general procedure for constructing the unstable 
fixed point and its scaling fields. In addition, the fact that we found overlapping domain 
of convergence in the particular computation presented here is an encouraging result but it 
does not guarantee that a similar overlap will be found for other models. 



II. THREE STEPS TOWARDS THE CALCULATION OF THE CRITICAL 

AMPLITUDES 

In this section, we consider a scalar model on a D-dimensional lattice with a lattice 
spacing oq. We assume that (3 < Pc and that the free energy density / is finite in the 
thermodynamic limit. We discuss the estimation of the critical amplitudes in the HT phase 
for x'^^\ the /-th derivative of — / with respect to an external magnetic field. For definiteness. 



we generalize the parametrization of Eq. (|1 . 1|) to 



X(') = - + - Pt''' + ...). (2.1) 

We assume the existence of a discrete RG transformation R, which can be performed in the 
following way. We first integrate the fields in blocks of side ha^ while keeping the sum of the 
fields in the block constant. We then rescale the sums of the fields by a factor 5(~2--D+r;)/2_ 
For the HM, h = 2^/^ and = 0. We assume that this RG transformation has one non-trivial 
unstable fixed point and an attractive HT fixed point. We first introduce general notations 
for the scaling fields and then discuss the three steps. 



A. Construction of the scaling fields 

A preliminary requisite is the construction of the two sets of scaling fields. We construct 
the scaling fields near the unstable fixed point, denoted y(d), in terms of the coordinates 
d in the directions of the eigenvectors of the linearized RG map at that fixed point. The 
boldface notations mean that the quantity is a vector. In the following we will consider finite 
dimensional approximations for such vectors. For small values of d, we have y(d) ~ d. One 
could call the d, the "linear scaling fields" or the "real fields" if they have a simple 
physical interpretation. If we denote the RG transformation in the d coordinates as R(d), 
and Xj as the eigenvalue in the j-th direction, we have by definition of the scaling fields 

y,(R(d)) = \,y,id). (2.2) 

In the following, we assume that the spectrum is real, positive, non-degenerate and that 
only Ai > 1. 
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Similarly, the HT scaling fields denoted y(h), can be constructed as expansions in the 
coordinates h in the directions of the eigenvectors of the linearized RG map at the HT fixed 
point. With notation analog to Eq. ( p.2| ), we have 

y,(R(h)) = A>,(h). (2.3) 

The d's can be expressed in terms of the h's by a shift followed by a linear transformation. 
In the following Un^i or ?/i(d„) will denote the value of the /-th scaling field after n iterations, 
and similar notations will be used for the HT variables. 



B. step 1 

The first step consists in expressing the x^^^ terms of the scaling fields of the HT fixed 
point y. In the vicinity of the HT fixed point, the effective lattice spacing becomes larger 
than the physical correlation lengths. This fixed point is attractive and as the number of 
iterations n becomes large, one can treat the blocks as almost isolated systems with a large 
volume 6"^. In the following we call x^n average value of the /-th power of the total 
spin in this volume minus its disconnected parts, divided by the volume. Assuming that 
lim„^ooXn'' = X*'''' is finite for (3 < Pc, we conclude that the subtracted average value of the 
l-th. power of the rescaled sum of all the spins scales like 5~({2+^-»?){'/2)-£')n. ^qj. ^ large and 
should be expressible as products of the y with product of eigenvalues 5~((2+-D-t;)(«/2)-d)_ 
Indeed, in explicit calculations for the HM model, we found that simple linear relations 
hold, namely, when n becomes large 

xi^^)ocy,(h„)(A;)-", (2.4) 

and consequently, 

^(2,) ^ K^i)y^{h,^) , (2.5) 

with denoting the initial values and K^'^^ a constant depending on the choice the scales 
of the coordinates h and easily calculable. This is discussed explicitly in subsection |V C 



for the HM. For the simplicity of the exposition, we will assume that this is also true for 
other models. Nevertheless it is straightforward to generalize the construction for the case 
of products of HT scaling fields. 



C. Step 2 

The expansion of ^g(hm) is expected to converge for hj„ small enough, however it might 
not be very useful or even meaningful near the unstable fixed point or the Gaussian fixed 
point. If 7 denotes the leading critical exponent for x^'^''\ it follows from arguments based 
on the linear RG transformation that 7*^^'^'' = - InAg/lnAi. Consequently, yj jjq is RG- 
invariant. We can then factor out the susceptibility into a singular part and a RG invariant 
part: 



^(2.)=^(.)[^^(d(h„,)))] 



[yi(d(h,„))]^'^^'y,(h. 



(2.6) 
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The factor in large brackets is RG-invariant and does not need necessarily to be evaluated 
for initial values of h. We can use the n-th iterate of these values h„ = R"(hj„) for any 
n and get the same answer. In particular, we can choose n is such a way that the flow 
is "in between" the two fixed points considered here, in a crossover region where the two 
expansions have a chance to be valid. One of the main results presented in this article are 
numerical evidences that this procedure actually works. In other words, that there exists 
an overlap between the domains where approximate expansions of the scaling fields scale as 
they should. 

The scaling properties of expansions are related to convergence issues which need to 
be discussed in the complexification of the construction. If an eigenvalue appearing in the 



defining equation for the scaling fields Eqs. (p.2|) and (|2.3|) , is such that \Xi\ ^ 1, it is clear 
that the only values that can be taken by \yi\ at a fixed point are and oo. Consequently, 
a detailed analysis of the fixed points for the complexification of the RG transformation 
can put restrictions on the domain of convergence of the scaling fields. Later, the inverse 
functions d(y) or h(y) will also be used. Their radius of convergence can be restricted by 
the study of the extrema of the original function. In the case of the one dimensional model 
of Section |^ such a study can be conducted easily and confirms the results of the numerical 
analysis. 

Near criticality, for the values of n in the crossover discussed above, the values of the 
scaling fields associated with the irrelevant directions are usually very small (unless very large 
initial values have been chosen) and can be treated perturbatively. In first approximation, 
we can consistently set the initial values of the irrelevant scaling fields to zero since they 
are multiplicatively renormalized. In this approximation, we describe the flow along the 
unstable manifold. A local calculation involving a scaling field corresponding to the the 
relevant direction is provided in Ref. following a procedure developed in Ref. to 
prove Feigenbaum conjectures. 

It should also be noted that since the RG considered here is discrete (it is constructed 
by iterating R an integer number of times), RG-invariant does not mean independent of 



hj„. There is room for log-periodic corrections, which have been first noticed by Wilson [15 



discussed in general in Ref. and observed in the HT expansion of the HM in Refs. |T3|,|n[. 
These corrections are studied for the simplified model in Section In our numerical study 
of the HM of Section |T1|, these corrections are too small to be resolved and will be ignored. 



D. Step 3 

In the previous step, we traded a difficult problem (the estimation of the initial values 
of y) for two simpler problems: the estimation of the RG invariant and the estimation of 
the initial values of y. Up to now, we have treated these initial values as free variables and 
constructed functions of these variables which depended only on the RG transformation. We 
now need to incorporate the information related to the actual initial values. This calculation 
depends on the type of models considered. For instance, for a Ising model, one expects 

yi = yi;i(/5, Y,,2{Pc -/?)' + ..., (2.7) 

and 
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(2.8) 



for / > 1. The Yi^^ are constants that we will evaluate numerically for the HM. The non- 
leading terms are responsible for the analytical corrections and are usually difficult to resolve 
numerically. A more complete construction of the initial values in terms of the "bare pa- 
rameters" used in field theoretical perturbative calculations (with flows starting near the 
Gaussian flxed point) is in progress ||29| but will not be discussed here. 



III. DYSON'S HIERARCHICAL MODEL 



In this section, we review the basic facts about the RG transformation of the HM to 
be used in the rest of the paper and discuss various way to obtain flnite dimensional trun- 
cations. For more detail, the reader may consult Ref. pl| , |23| and other papers quoted in 
the introduction. The energy density of the HM has two parts. One part is non-local (the 
"kinetic term" ) and invariant under a RG transformation. Its explicit form can be found, for 
instance, in Ref. ||21| or in Section II of Ref. ||3^. The other part is a sum of local potentials 
given in terms of a unique function V{(j)). The exponential e"^'^'^-' will be called the local 
measure and denoted 14^0 (0)- For instance, for Landau- Ginsburg models, the measures are of 
the form WQ^cj)) = e~'^'^ , but we can also consider limiting cases such as a Ising measure 
VFo(0) = (^(0^ ~ !)■ Under a block spin transformation which integrates the spin variables in 
"boxes" with two sites, keeping their sum constant, the local measure transforms according 
to the integral formula 



+u 



+lg{/3/2)(c/4)"+V2 



Wr,. 



+ 



(3.1) 



where f3 is the inverse temperature (or the coefficient in front of the kinetic term) and Cn+i 
is a normalization factor to be fixed at our convenience. 
We use the Fourier transform 



dk 
2^ 



(3.2) 



and introduce a rescaling of /c by a factor m/s", where u and s are constants to be fixed at 
our convenience, by defining 



(3.3) 



In the following, we will use s = 2/^0 with c = 2^^^/^. This corresponds to the scaling of 
a massless gaussian field in D dimensions. Contrarily to what we have done in the past, we 
will here absorb the temperature in the measure by setting u = With these choices, 
the RG transformation reads 



Rn+i{k) = Cn+iexp 



1^ 
2dk^ 



Rn 



\fck^ 



(3.4) 
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We fix the normalization constant C„ so that -Rn(O) = 1. For an Ising measure, Ro{k) = 
cos{^/J3k), while in general, we have to numerically integrate to determine the coefficients 
of Ro{k) expanded in terms of k. 

If we Taylor expand about the origin, 

oo 

/?„(fc) = l + 5]a„,;F, (3.5) 



1=1 



the finite-volume susceptibility reads 



X^ = -2^(^) . (3.6) 



The susceptibility x is defined as x = linin^oo Xn ■ For P larger than Pc, the definition of x 
requires a subtraction (see e. g., Ref. |^0[ for a practical implementation). In the following, 
we will only consider the HT phase (/3 < /3c). The explicit form for a^+i,; in terms of a„_; 
reads 

O'n+1,1 — 5 1<J- ' J 

Unfi 

where 

_^ (-l)-(2(/ + .))! 

Un,l - 2^ „2(/+i)„-lf97^|l «n,pan,g • l<3-ej 

1=0 ^ p+q=l+i 

To study the susceptibility not too far from the HT fixed point, we can expand x 
terms of /5. Since we choose the scaling factor u so that P is eliminated from the recursion, 
we find that Oq,; oc pK From the form of the recursion, Eq. (|3.8|) , we can see that an,i will 
always have as the leading power in its HT expansion (since p + q > I)- We define the 
coefficients of the expansion of the infinite- volume susceptibility by 

oo 

X{P) = E bmP"" . (3.9) 

m=0 

We define = bm/bm^i, the ratio of two successive coefficients, and introduce quantities 
I9| , called the extrapolated ratio (Rm) and the extrapolated slope (Sm) which will be used 



later in a more general context and are defined by 

Rm = mrm - (m - l)rm-i , (3.10) 

and 

Sm = mSm - (m - l)Sm-l , (3.11) 

where 

g,,.^ (3.12) 

mrm - im- l)r^_i 
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is called the normalized slope. If we calculate Sm for the HM, we find oscillations illustrated 
in Refs. 



The HT expansion can be calculated to very high order, however, due to the amplifica- 
tion of some subleading corrections by the estimators, this is an inefficient way to obtain 



information about the critical behavior of the HM. In Ref. [|23|, it was found that one can 
obtain much better results by neglecting the contribution of the an,i when / > Imax-, with 
Imax much smaller than the order of the HT expansion. As an example, one can calculate 
the 1000-th HT coefficient of x with 16 digits of accuracy using only 35 terms in the sum. 
This is equivalent to consider the polynomial approximation 

^max 

R^{k) ~ ' (3-13) 

1=0 

for some integer Imax- There remains to decide if one should or not truncate to order /c^Uaa: 
after squaring This makes a difference since the exponential of the second derivative 
has terms with arbitrarily high order derivatives. Numerically, one gets better results at 
intermediate values of Imax by keeping all the terms in R^. In addition, for the calculations 
performed later, the intermediate truncation pads the "structure constants" of the maps 
(see sec. 0) with about fifty percent of zeroes. A closer look at Section 0, may convince 
the reader that not truncating after squaring is more natural because we obtain correct (in 
the sense that they keep their value when Imax is increased) structure constants in place of 
these zeroes. We have thus followed the second possibility where we truncate only once at 
the end of the calculation. With this choice 

"(-im[+^ ^ 

Compared to the HT expansion, the initial truncation to order Imax is accurate up to 
order p'-"^^^. After one iteration, we will miss terms of order ^''"'»^+^ but we will also generate 
some of the contributions of order z^^'"""" (but not all of them). After n iterations we generate 
some of the terms of order z^^"'"""" as in super convergent expansions. Rigorous justifications 



of the polynomial truncation can be found in Ref. p2| 



IV. A ONE- VARIABLE MODEL 

Before attacking the multivariable expansions of the scaling fields, we would like to 
illustrate the main ideas and study the convergence of series with a simple one variable 
example which retains the important features: a critical temperature, RG fiows going from 
an unstable fixed point to a stable one, and log-periodic oscillations in the susceptibility. 

In order to obtain a simple one- variable model, we first consider the Imax = 1 truncation 
ofEq. (|3l4l ): 



(c/2)a„,i - (3cV8)a^,i 
""-^^'^ 1 - (c/2)a„,i + (3cVl6)a2 " ^^''^ 



Expanding the denominator up to order 2 in a^+i i and using Eq. (|3.6| ), we obtain 
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Xn+l = Xn + j[£)^ Xl- (4.2) 



This approximate equation was successfully used in Ref. p3| to model the finite-size effects 
and was used as the starting point for a study of the scaling field in Ref. . If we expand 
X (the limit of Xn when n becomes infinite) in /?, and define the extrapolated slope, Sm, as 
in Eq. (|3.11| ), we see oscillations in Fig. ^ quite similar to those in the HM [13,14|. Using a 
rescaling discussed in Ref. [|18| and the notation ^ = c/2, the map can be put in the canonical 
form 

hn+l=^hn + {l-Ohl. (4.3) 

We call this map the "/i-map" . It has a stable fixed point at and an unstable fixed point 
at 1. We recover the susceptibility as: 

Xn = (4.4) 

We can expand the map about the unstable fixed point, hn = 1- Using the new coordinate 
(i„ = 1 — hn and the notation A = 2 — ^, we obtain the "dual" map 

rf„+i = Xdn + (1 - \)dl , (4.5) 

with the starting value do = 1 — (3/ (3c- We call this map the "rf-map" and we can think 
of d as being the distance to the critical point. We can construct a function d such that 
d{iin) = dn by plugging an expansion in into the equation 

d{\yn) = \d{yn) + (1 - \)d\yn) . (4.6) 

Similarly, we can construct the inverse series and obtain d{yn)- This allows us construct dn 
in terms of do'. 

dn = y-\XMdo)) . (4.7) 

Because of the duality between our two maps, we can easily reproduce all of the above results 
of the (i-map for the /i-map and express hn in terms of a HT scaling field ijn- 

We now turn to the three steps. Note that step 3 is not necessary: the knowledge of 
y{do) provides the initial value of y as a function of /3. Step 1 is straightforward. As shown 
in Ref. |T^, in the infinite n limit, 

X- = f = ^. (4.8) 

ho ho 

The only difficult part is step 2. The susceptibility can be written as 

^ = (l-do)(y(rfo))^ ■ ^^'^^ 

with the RG invariant G = yoVo- Due to the discrete nature of our RG transformation, 
is not exactly a constant. If expressed in terms of ln(?/o), © is a periodic function of period 
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InA. However, for A not too close to 2, the non-zero Fourier modes are very small. Note also 
that the apparent singularity when (ig — > 1 is exactly canceled by y{do)'^ by virtue of Eq. 
( [4.11| ) discussed below. We now give empirical results concerning the large order behavior 
of of the expansions of y{d), y{h) and their inverses. 

We first consider d{y) = y + Z]j^2 ^iv''- For all tested values of 1 < A < 2, we obtained 
very good linear fits of ln|si_i/s;| versus ln(Z), for / large enough. Thus for large /, the 
coefficients obey the approximate rule 



Sl+l 



CI'' , (4.10) 



where we find that C is always of order 1 (in fact, 0.9 < C < 1) and < < 1. Using 
iteratively this formula, we find that the coefficients decrease like C^'(Z!)~'^ and consequently 
the d{y) should be an entire function. The numerical values of k are given in Fig. ^. 

We then consider h{y) and again examine the ratios of successive coefficients, 
We find the ratios fiatten to constant values, for large enough /, indicating a finite radius of 
convergence. The radius get smaller and vanish as ^ approaches zero as shown in Fig. ^. 
Note that for any value of ^ tried, we found good evidence that lirrin^aoSn = —1-5, indicating 
a (y — ycY^"^ behavior. This is consistent with the existence of a quadratic minimum for the 
inverse function discussed below. 

We now turn to the inverse functions starting with y{d). As (i — > 1, we reach the HT 
fixed point and we expect the convergence of the series to break down in this limit. We find 
empirically from the analysis of ratios that for all 1 < A < 2, |/((i) converges in the region 
< (i < 1. The analysis of the extrapolated slope for various A gives convincing evidence 
that the main singularity has the form 

y(rf) ~ (1 - . (4.11) 

This can be seen with short series when A is close to one and requires larger and larger series 
as A gets close to 2. Illustration of these properties for a particular value of A are shown in 
Figs. ^ and ^ where the estimators are compared with those of the HM. 

Finally, we discuss y{h) which can be seen as a high-temperature expansion. We found 
clear evidence that the ratio of coefficients tj/ti+i approaches 1 for large I. For smaller values 
of ^, it takes larger order to reach this limit. A detailed study shows that if we continue 
y{h) for negative values of h using the series expansion, the function develops a quadratic 
minimum at some negative value of h. The absolute value of y at that value of in all 
examples studied, reproduces accurately the radius of convergence of the inverse function. 
As one may suspect by looking at Fig. ||, the analysis of the extrapolated slope is intricated. 
However, if we calculate enough terms and if A is not too close to 2, we get approximate 
results which are consistent with a main singularity of the form 

y{h)-{^-h)-^ ■ (4.12) 
For instance, just by looking at the asymptotic behavior of Fig. 0, one can see that 7 ~ 



1.4677, as expected, with errors of the order 10 ^. It is interesting to note the duality |1S 
between Eqs. (|TTp and ( CT )- 

All the results concerning the convergence and the singularities have a simple interpreta- 
tion. The finite radius of convergence of y and y is due to the other fixed point which cannot 
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be located inside the domain of convergence. In this simple example there are exactly two 
fixed points in the complexification of the map and this concludes the discussion. On the 
other hand, the finite radius of h is due to a minimum of y at negative values of h while 
such a minimum does not appear for d. 

As we have seen above, ynUl is independent of n. We can thus pick n such that we are 
just in the crossover region and both expansions have a reasonable chance to be accurate. 
In order to test the accuracy of the approximations yapp{d) (series expansion up to a certain 
order) for various n, we have prepared an empirical sequences of starting with do = 10~*. 
We have then tested the scaling properties by calculating 

i^n=|[2/appK)/(2/0A")]-l| , (4.13) 

where yo was calculated with 16 digits of accuracy by using enough terms in the expansion 
of y{do). For double precision calculations, optimal approximations are those for which 
Dn — 10~^^. For such approximation, the scaling is as good as it can possibly be given 
the accuracy of yo. Indeed, due to the peculiar way numerical errors propagate ||31|, one 
does not reach exactly the expected level 10^^^. We can define a similar dual quantity by 
replacing d hj h and y by y. In this case, yo is estimated with the same accuracy as yo by 
stabilizing y{hn)/^^, for large enough n. 

We have performed this calculation for A = 1.1, 1.5 and 1.9. The conclusions in the 
three cases are identical. For n large enough, the Dn of y starts increasing from 10~^^ until 
it saturates around 1. By increasing the number of terms in the expansion, we can increase 
the value of n for which we start losing accuracy. Similarly, for n low enough, the Dn of y 
starts increasing etc... We want to know if it is possible to calculate enough terms in each 
expansion to get scaling with some desired accuracy for both functions. The answer to this 
question is affirmative according to Fig. |^ for A = 1.5. One sees, for instance, that with 
10 terms in each series, we have scaling with about 1 part in 1000 near n = 45 for both 
expansions. The situation can be improved. For 70-70 expansions, an optimal accuracy is 
reached from n = 44 to 46. For the other values of A quoted above, similar conclusions are 
reached, the only difference being the optimal values of n. 

Another evidence for overlapping convergence is that we can stabilize the RG invariant 
6 for a certain range of To evaluate 0, we use the series expansions for y and d, cutting 
each off at some order: 

m m 

y(i-%))^E^^(i-E^.vr, (4.14) 
i=i j=i 

where si and ti are the /th coefficients in the d and y series, respectively. We have found 
that, given a fixed value of m + rh, the most accurate values for are obtained when m ~ m. 
In Fig. ^, we show calculated by keeping 50 terms each in the expansions for y and y. 
The result is plotted against ln(?/). We used ^ = 0.1, which makes the oscillations much 
larger than, for example, near to ^ = 2~^/^. Near the fixed points, we need more terms in 
the appropriate series to get accurate results. 

We can study the oscillation we see in by looking at its Fourier expansion. Since is 
periodic in Inyo, we can express 

eiyo) = yM^ - d{yo)) = E ape'^'-'-y , (4.15) 
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where uo = 27r/lnA. The coefficients are given as 



Op = ^'"'"'"^(1 - d^y))dy > (4-16) 

As an example, we calculated Oq for ,^ = 0.1 where the oscillations are not too small. The 
choice of the interval of integration can be inferred from Fig. ^. If we had infinite series, the 
function would be exactly periodic. For finite series, we see that Ua cannot be too large or 
too small. For intermediate values, we obtain cq — 6.06676. Proceeding similarly, we were 
able to resolve the next two Fourier modes. For reference, the magnitude of 02 is about 
2 X lO"^'^. Using B ~ ao together with Eq. (|4.9| ), we obtain the leading amplitude together 
with the analytical corrections coming from the nonlinear tems in y{d). 



V. SCALING FIELDS IN THE HIERARCHICAL MODEL 

A. Construction of the scaling fields 

For notational convenience, we first rewrite the unnormalized recursion given in Eq. 
( p.l4| ), using the "structure constants": 



with 



[ , otherwise . 

These zeroes can be understood as a "selection rule" associated with the fact that an,i is 
of order (3^ as explained in Section If we follow the truncation procedure explained in 
Section |In| , the indices simply run over a finite number of values. We use "relativistic" 
notations. Repeated indices mean summation. The greek indices indices /x and v go from 
to /max, while latin indices i, j go from 1 to Imax- Obviously, F^'' is symmetric in fi and z/. 
By construction, a„ = 1 and we can write the normalized recursion in the form: 

- ^ —;rj , lo-ij 

with A^; = 2ro^ . 

We then expand the basic map of Eq. ( p.3|) about the two fixed points of interest, choosing 
coordinates such that the matrix associated with the linearized RG transformation becomes 
diagonal. This matrix is not symmetric and the relations of orthogonality and completeness 
are left invariant under rescaling of any right eigenvector by a nonzero constant together 
with a rescaling of the corresponding left eigenvector by the inverse. In the following, we will 
fix this ambiguity by requiring, in analogy with Section |^ that the "other" fixed point (i.e. 
the one not located at the origin by construction) be at (1, 1, . . .) in the new coordinates. 

The HT fixed point is at the origin of the coordinates in Eq. ( p.3| ) and all we need 
to do is diagonalizing the linear form This is quite simple because it is of the upper 



14 



triangular form. The eigenvalues are just the diagonal terms. From Eq. ( ^.2|) , one sees that 
Ith diagonal term is given as 



A, = 2(c/4)" . (5.4) 



This spectrum was obtained in Lemma 3.3 of Ref. ||2T[ with a different method. This 
spectrum has a simple interpretation which will be discussed below. The upper diagonal 
form implies that the / — th right eigenvector has only its / first entries non-zero. This means 
that if we truncate to ai^^^, we are simply truncating to a subspace spanned by the first Imax 
eigenvectors. This is an interesting reinterpretation of the original polynomial truncation 
which can be applied for other models. Introducing new coordinates hi, so that a„,i = ipihn,r 
with the matrix of right eigenvectors, the RG transformation has the form 

\ h -L /\P1h h 

/\j.ibn,r ~r ij. ' i'n,p' <'n,q /- -\ 

iln+l,r = r—, '■, ■ l^-^J 

1 + AP/i„,p + Ag^/i„>„,, 

Note that the form of the eigenvectors guarantees that hn^i is of order (5^ . This can be seen 
by inverting the linear change of variable using the matrix of left eigenvectors. Due to the 
upper-diagonal form of the linear transformation, the second left eigenvector has its first 
entry equal to zero, the third its first two entries etc... . 

Near the nontrivial fixed point, we first use accurate values of the fixed point to 
bring the fixed point at the origin. The eigenvectors are then calculated numerically using 
truncated forms of the linearized RG transformation. There is no exact closure as in the HT 
case, however the first eigenvalues stabilize rapidly when Imax increases. There is only one 
eigenvalue larger than 1. For instance the numerical values for c = 2^/^ are Ai = 1.42717 . . ., 



A2 = 0.85941 .... A more complete list is given in Ref. In summary, we can choose a 
system of coordinate di where the unstable fixed point will be at the origin of the coordinate 
and the HT fixed point at (1, 1, . . .), and such that the RG transformation has the form 

\ d -I- AP'^r/ d 

1 -|- A^dn.p + /S.^'^dn^pdn^q 

We can express canonical coordinates (linear scaling fields) in terms of the nonlinear 
scaling fields: 

dn,r = ^ ^r,iii2...yn,iyn,2 ■ ■ ■ i (5-'^) 
11,12, ■■■ 

where the sums over the i's run from to infinity in each variable and yn,i = ^^yo,i- Using 
the notation i = (ii, 12 ■ ■ ■) and the product symbol, we may rewrite the expansion as 



m 



Using the transformation law for the scaling fields, we have 

dn+l,r = J2 ^r,i HiXmymY"' ■ (5.9) 
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Each constant term, 5^,0,0,..., is zero, as the scahng fields vanish at the fixed point. From 
Eq. (|5.6| ), we see that all but one of the linear terms are zero for each value of r. The 
remaining term is the one proportional to the rth scaling variable. We take these coefficients 
to be 1, so that the dn^r — yn,r for small ijn^r- For the higher-order terms, we obtain the 
recursion 

X]j+k=i('^r'''^P,j'^<?,k -^rj Tim ^In ■'^p,k) X]j+k+l=i ^r,} Tim ^rn ^0 "^Pik-^g,! TC[\ 

(n.A--A.) • ^'-'^^ 

The calculation can be organized in such way that the r. h. s. of the equations are already 
known. This will be the case for instance if we proceed order by order in iq, the degree of 
non-linearity. This expansion may have small denominator problems. However, as discussed 
in the introduction, using numerical values of the eigenvalues as calculated in [^, we did 
not find spectacular cancellations between the two terms entering the denominator. 

We can likewise expand each hn,r in terms of scaling fields yn,i,yn,2, ■ ■ ■■ The derived 
recursions are identical in form to those derived above. From Eq. (^.4|), one sees that the 



denominator will vanish for some equations. This question is discussed in p6[ where it is 
shown that to all order relevant for the following calculation, a zero denominator always 
comes with a zero numerator. 

One can likewise find expansions of the scaling fields in terms of the canonical coordinates, 
by setting 

i m 

and requiring that when d„ is replaced by d„_|_i, the function is multiplied by A^. Since d„+i 
has a denominator, it needs to be expanded for instance in increasing order of non-linearity. 
A simple reasoning shows that exactly the same small denominators as in Eq. ( [5.10[ ) will be 
present in these calculations. The same considerations apphes for y. 



B. Practical implementation 

We have calculated an empirical series of a„j with c = 2^/^, an initial Ising measure and 
(3 = (3c — 10~* . Detail relevant for this calculations can be found in Refs. [^3,24,3^. The 
calculations have been performed with Imax = 30, a value for which at the (3 considered, the 
errors due to the truncation are of the same order as those due to the numerical errors. These 
errors are small enough to allow a determination of the susceptibility with seven significant 
digits if we use double precision. 

We now discuss the fiow chronologically. Our choice of (3 (close to (3c) means that we 
start near the stable manifold. After about 25 iterations, we start approaching the unstable 
fixed point and the linear behavior dn+i,i — \idn becomes a good approximation. During 
the next 20 iterations, the irrelevant variables die off at the linear rate and at the same 
time the flow moves away from the flxed point along the unstable direction, also at the 
linear rate. At = 47 we are in good approximation on the unstable manifold and dn,2 
becomes proportional to d^i- In other words, the non- linear terms are taking over. At this 
point, we can approximate the dn,i as functions of yi only: dn,i — (ii(A5'?/i, 0, 0, . . .). This 
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approximation is consistent in the sense that if ?/2 = at n = then it is also the case for 
all positive n. 

We have calculated diiiji, 0, 0, . . .) up to order 80 in yi using Eq. ( p.lO| ). There cannot 
be any small denominators in this restricted case. We have then inverted di, now a function 
of Hi only, to obtain yi{di). Given the empirical dnA we then calculated the approximate 
yn,i and then used the other functions diiyi) (with / > 2) calculated before to "predict" 
dn^i- Comparison with the actual numbers were good in a restricted range. For n = 49, the 
relative errors were less than a percent. They kept decreasing to less than one part in 10,000 
for n = 54 and then increased again. It will be shown later that this corresponds to the fact 
that when yi becomes too large (a value of approximately 3.7 first exceeded at n = 57), the 
series expansion of di seem to diverge, unlike the one-variable model for which d{y) is an 
entire function. 

The quality of the approximation between n = 45 and n = 55 can be improved by 
treating y2 perturbatively. We have expanded 

di{y,,y2, 0, . . .) ~ rff ^2/1, 0, 0, . . .) + y2di'\yu 0, 0, . . .) , (5.12) 

with d\^^ up to order 30 in yi. This allows us to obtain the first order expression: 

UHd{yl,y2,0,...)))-Giy^)+y2H{y,) . (5.13) 

Note that expansions at both fixed points are involved (one for yi and one for d) in this 
equation. When finite series are used, the approximation will only be valid in the crossover 
region. Near n = 55, the presence of the HT fixed point starts dominating the flow but 
we are still far away from the linear regime. We have taken these non-linear effects in 
^i(h) into account by calculating it up to order 11 in f3. This is a multi- variable expansion. 
Recalling the discussion about the HT expansion in Section |in| and the properties of the 
eigenvectors of the linearized RG transformation about the HT fixed point, we can count 
the number of terms at each order in (3. At order two, we have hf and h2, but since the 
linear transformation is diagonalized, h2 will only appear in y2{h.) with coefficient 1. It is 
easy to see that at order /, one needs to determine p{l) — 1 coeffificients, where p{l) is the 
number of partitions of /. It has been known from the work of Hardy and Ramanujan that 



exp(7r./2//3) 

It seems thus difficult to get very high order in this expansion. In order to fix the ideas, 
there are 41 terms at order 11, 489 at order 20 and 13,848,649 at order 80. 

As we will explain below, the expansion up to order 11 has a sufficient accuracy to be 
used starting at = 55. It also provides optimal (given our use of double precision) results 
for 77, > 60. As n increases beyond 60, one can see the effect of each order disappear one after 
the other as discussed in |VD . Finally, the linear behavior becomes optimal near n = 130. 



This concludes our chronological discussion. 

C. Step 1 

From Eqs. (|3.6|) and remembering that we have absorbed (3 in the a„ , we obtain 



17 



Xn = -{2/(3WM'2/cr , (5.15) 

where if) is the matrix of right eigenvectors. For n large enough, the hnear behavior apphes 
and the hn^r get multiphed by 2(c/4)^ at each iteration. In the large n limit, only the r = 1 
term survives and consequently, 

X = -(2//5)^Jlim™/in,i(2/c)" . (5.16) 

Using the same method as in the one- variable model, we can in the limit replace hn^i by y„ i 
and obtain 

X = -{'^/m\yo,i ■ (5.17) 

For reference, in the case c = 2^^'^ and with the normalization of the eigenvectors discussed 
above, ~ —0.564. Also note that since for all /, < < 1, all other monomials in the yi 
go to zero faster than yi. 

One can calculate the subtracted 2g-point function following the same procedure. As 
shown in [^,0 , they can be expressed in terms of aq^n and higher powers of the an,i with / < 



q. Following the procedure described above, these quantities can then reexpressed in terms 
of y. We need to identify the leading term in this expansion. By rescaling (— l)''(2g)!aq „ 
by (4/c)" we obtain the average value of the 2g-th power of the main spin (sum of the 2"' 
spin variables (px)- In the symmetric phase, this quantity scales like 2''"' when n increases. 
However, the subtracted version of this quantity (which is generated by the free energy) is 
expected to scale like 2^. In other words, if we assume that the free energy density is finite 

ttq^n - (subtractions) oc [2(c/4)'?]" . (5.18) 

One clearly recognizes the spectrum of Eq. ( |5.4|) and expects that the leading term is 

o-g.n — (subtractions) oc yn^q ■ (5.19) 

In order to prove this conjecture by direct algebraic methods, one needs to show that the 
the non-linear terms which scale faster than ijn^q are canceled by the subtraction. For 
instance for the subtracted four-point function and c = 2^/^, > A2 and X\ > A2. We 
have checked that the corresponding terms {yiY and (yi)^ disappear with the subtraction. 
We have conducted similar checks for the 6 and 8 point functions p9| and found similar 



cancellations. It should be noted that a rigorous proof that the mechanism works in general, 
would imply the finiteness of the free energy density and hyperscaling (as defined in [pO[]). 
The practical calculation of the subtracted quantities is made difficult by the fact that as 
n increases, the "signal" becomes much smaller than the "background" (the unsubtracted 
part). This requires the use of adjustable precision arithmetic. In the following, we will only 
discuss the 2 point function (susceptibility). 

D. Step 2 

We rewrite the susceptibility as 

X^(1.127853//5)e(yo,l)'^ (5.20) 
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with the RG invariant 9 = yo,il/o,i = yn,iy1,i- We first constructed yi{di) by neglecting the 
effects of the irrelevant directions, as explained in subsection VB| . In order to provide a 
comparison, we have calculated the same y{d) for the one-variable model with A = 1.2573. 
In the following, we call this model the "simplified model" (SM). For this special value 
of A, the critical exponents 7 of the two models coincide with five significant digits. We 
have good evidence that both series have a radius of convergence 1 as indicated by the 
extrapolated ratio defined in Eq. ( ^.l(j| ) reaching 1 at an expected rate (Fig. Similarly, 
their extrapolated slopes seem to converge to the same value I/7 — 1 ~ —0.23026. . . as 
shown in Fig. |^. In conclusion, the function y{d) for the SM is a reasonably good model to 
guess the asymptotic behavior of yi{di). Remembering the explanations of Section |rV|, this 
indicates the absence of other fixed points in the vicinity of the unstable fixed point. 

For the inverse function di{yi), the situations more complicated as shown in Fig. |^. 
The quantity plotted in this figure will be used to discriminate between a finite and an 
infinite radius of convergence. If 1 6m I ~ R~"^ as for a radius of convergence R, then we have 
ln{\bm\)/fn ~ —ln{R) + A/m for some constant A. On the other hand, if \bm\ ~ (m!)~" 
as for an infinite radius of convergence, then we have ln{\bm\)/iTL ~ — Q;(ln(m) — 1). In the 
following, we will compare fits of the form Ai + A2/m and i?iln(m) + B2. We first consider 
the case of SM, where according to the our study in Section [rV| , expect an infinite radius of 
convergence. This possibility is highly favored as shown in Fig. |10[ One sees clearly that 
the solid line is a much better fit. The chi-square for the solid line fit is 200 times smaller. 
In addition Bi ~ —B2 as expected. 

The analysis for the HM is more delicate. One observes periodic "dips" in Fig. ^ 
which make the ratio analysis almost impossible. We have thus only considered, the "upper 
envelope" by removing the dips from Fig. The fit represents an upper bound rather than 
the actual coefficients. The fits of the upper envelope are shown in Fig. |TT]. The possibility 
of a finite radius of convergence is slightly favored, the chi-square being 0.4 of the one for 
the other possibility. Also, the second fit does not have the Bi ^ —B2 property. From 
Ai ~ —1.32, we estimate that the radius of convergence is about 3.7. This means that if we 
want to write some analytical formula for the fiows by first calculating the initial values of 
the scaling fields and then calculating d at successive iterations by using their expression in 
terms of the scaling fields, we will have to switch variables in the crossover region. 

As explained above, the approach of the HT fixed point is intrinsically a multivariable 
problem. For this reason, the calculation of D^, defined in Eq. ( |4.13|) , that tests the quality 
of scaling, will be our main tool of analysis. In the following, we limit the discussion to 
^i(h) which enters in the susceptibility. We have calculated yi(h) in terms of 11 variables, 
up to order 11 in 13. As in section |V], we will use an empirical series a„ i, calculate the 
corresponding hn,i an plug them in the scaling fields and check the scaling properties. This 
empirical series was calculated with an initial Ising measure and (3 = (3c — 10~^ (see Ref. 



Again we define a quantity as in Eq. ( [4.13[ ) which is very small when we have 
good scaling and increases when the approximation breaks down. The results are shown in 
Fig. for successive orders in (3. The solid line on the right is the linear approximation. It 
becomes optimal near n = 130. The next line (dashes) is the second order in f3 expansion. 
It becomes optimal near n = 90. Each next order gets closer and closer to be optimal near 
n = 60. The last curve on the left is the order 9 approximation. It is hard to resolve the 
next two approximations on this graph. 
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The asymptotic value is stabilized with 16 digits and one may wonder why we get only 



scaling with 14 or 15 digits in Fig. |12|. The reason is that we use empirical data and that 



numerical errors can add coherently as explained in Ref. |^T[|. This can be seen directly by 
considering the difference between two successive values of Dn. A detailed analysis shows 
that the numerical errors at each step tend to be negative more often than positive, and 
consequently there is a small "drift" which affects the last digits. 

We can now look at the Dn defined as in Eq. ( [4.13| ) for yi and iji together in Fig. |TB. 



If we use an expansion of order 5 in /5 for y and of order 10 in cii for yi, we can get scaling 
within a few percent for both variables at n = 54. We can go below 1 part in 1000, with an 
expansion of order 11 in /3 and order 80 in di. At this point, the main problem is that the 
effects of the subleading correction makes the scaling properties worse when n < 57 and n 
decreases. One can improve the scaling properties by taking the effects of y2 into account. 
A detailed study shows that one can estimate the subleading effects between n = 40 and 
n = 45 as 

~ 7.2778 X 10"*^ + 3.2 x 10"^ x (5.21) 

It is thus possible to get a function scaling better by subtracting these correction. This 
improve the scaling properties by almost one order of magnitude near n = 54 and by almost 
two order of magnitude near n = 45. From Fig. |T^, we see that the combined scaling is 
optimal near = 54 which corresponds to an approximate value of 2 for yi. 

Using Eq. ( p. 121 ), we can take into account the first order correction in ?/2- After simple 
manipulations, we can rewrite the RG invariant as 

e-Ci + C2y2,o{yifi)^ , (5.22) 



with Ci = G{yi)yi and H{yi)yJ . The functions G and H are defined by Eq. ( p. 131 ). 



They rely on both expansions used and consequently they are only valid in a crossover region. 
Using explicit forms of Ci and C2 as a function of yi, we observe a plateau for each function 



which are shown in Fig. Using the flattest part of the plateau to estimate the constant 
we obtain Ci ^ 1.46416 and C2 — 1.663. Note that these two numbers are dependent 
of the choice of the scales for the scaling fields, but independent of the initial conditions. 
Consequently, if everybody agreed on the scales, these quantities could be called universal. 



E. step 3 

There remains to estimate the initial values of the yi and 1/2 • This step will done nu- 
merically from empirical values of d„. First, we obtain a rough estimate of ?/2,n from (i2,n at 
values of n where the linear approximation is good. Plugging this value in Eq. ( |5.12| ) for 



/ = 1, and inverting to get yn,i as a function of dn,i- Dividing by A" we get an estimated 
value of 1/0,1 with a plateau of about 10 iterations where the value is stabilized with 6 digits. 
Repeating for various temperatures we were able to determine the leading coefficient in Eq. 
(p.7|): Yi-i ~ 0.72782. Using Eq. ( |5.12| ) but for d2, together with our previous estimates of 
yi-n, we obtain ¥2-0 ~ —0.565. Subleading coefficients are difficult to resolve because 2A ~ 1 
and the nonlinear contributions in 7/2 are of the same order as the analytical corrections. 
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This conclude our approximate treatment of step 3. Using Eqs. ( |5.20| ) and ( ^.22|) we 
obtain the usual parametrization of the susceptibility of Eq. (LI) with Aq ~ 2.1162 and Ai ~ 
—1.196 in very good agreement with a fully numerical determination of these amplitudes 
0. 



VI. CONCLUSIONS 

In summary, we have shown with two examples that the 3 steps advocated in Section 
m lead to an accurate determination of the leading and subleading critical amplitudes. We 
provided numerical evidence that the formal expansions of the scaling fields associated with 
the HT and unstable fixed point have reasonable convergence properties and scale properly in 
overlapping domains. The determination of the initial values of the scaling fields associated 
with the unstable fixed point (third step) have been obtained numerically in the case of an 
initial Ising measure. Putting everything together, we were able to confirm numerical results 
obtained earlier. 

Analytical methods are now being used E^ to consider initial measures of the Landau- 



Ginzburg type in the vicinity of the Gaussian fixed point. For such initial measures, we can 
use perturbation theory in the quartic (or higher orders) coupling constant to construct the 
scaling fields associated with the Gaussian fixed point. Interestingly, one could use some 
of the methods developed here to interpolate between the Gaussian fixed point and the 
unstable fixed point. The completion of this task will allow us to give analytical expres- 
sions for the renormalized quantities in terms of the bare ones, which is the notoriously 
difficult problem that has to be faced by a field theorist. We are planning to extend the 
calculations performed here, to higher order derivatives of the free energy with a non-zero 



magnetic field and check explicitly amplitudes relations appearing in the literature p^-p4 



Another interesting question that could be addressed within this context is the crossover 



from classical to critical behavior |]7|,|9|,pl)|| . The completion of these projects will provide a 
detailed comparison between general RG expectations and their practical realization for the 
HM. 

The hierarchical approximation used in this article has allowed us to calculate large 
order expansions for the scaling fields. The fact that the general ideas advocated have 
worked properly means that one should now attempt to apply them to nearest neighbor 
models where similar calculations would be more time consuming. The examples we have 
in mind are spin models in three dimensions and asymptotically free theory such as the 
0{N) spin models in two dimensions or lattice gauge theory in four dimensions. In all cases, 
there exists some advanced technology for the weak and strong coupling expansions but the 
question of the interpolation has only been studied with the Monte Carlo method 

We expect that some of the simple features found in the study of the HT fixed point of 
the HM will generalize to nearest neighbor models. First, the fact that the HT scaling fields 
are linearly related to the successive derivatives of the free energy. Second, the fact that the 
restriction to a finite number of eigenvectors of the linearized RG transformation near the HT 
fixed point can be used to obtain improved HT expansions such as the polynomial truncation 
used above. However, the most difficult task remains a construction of the unstable fixed 
point and the RG fiows in its vicinity, with a control comparable to the case the hierarchical 
model. This is a challenging problem for the future. 
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FIGURES 




RG Flow 



G:Gaussian FP 
W: Wilson's FP 
HT: HT FP 



FIG. 1. Schematic representation of a RG flow starting near the Gaussian flxed point, passing 
near Wilson's fixed point and ending at the stable high-temperature fixed point. The circles 
represent the domains of validity of expansions of the scaling fields near the three fixed points. 
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FIG. 2. The extrapolated slope (S'„ 
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FIG. 3. The exponent k defined in Eq. ( [4.10| ) as a function of A. 
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FIG. 4. Radius of convergence of h{y) as a function of A = 2 — ^. 
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FIG. 5. Departure from scaling £)„ defined in the text, for y (curves reaching 1 to the right) 
and y (curves reaching 1 to the left). In each cases, we have considered approximations of order 
10 (dot-dashes), 30 (dots), 50 (dashes) and 70 (solid line). The value of A is 1.5. 
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FIG. 6. The invariant function O, calculated at ^ = 0.1, and plotted against In(yo)- 



30 




I I 1 I 1 I 1 I 1 I 1 ILl 1 I 1 I 1 I 1 

10 20 30 40 50 60 70 80 90 

Order m 



FIG. 7. The absolute value of the difference between the extrapolated ratio and 1 for the HM 
(empty boxes) and the SM (full circles), as a function of the order. 
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FIG. 8. The extrapolated slope Sm for the HM (empty boxes) and the SM (full circles) as a 
function of the order . 
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FIG. 9. Logarithm of the absolute value of the coefficients of the expansion of cZi(yi) divided 
by the order, for the HM (empty boxes) and the SM (circles). 
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FIG. 10. Comparison of fits of the form Ai + A2/m (dots)and i?iln(m) + B2 (solid line) with 
the data provided in Fig. ^ for the SM (circles). 
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FIG. 11. Comparison of fits of the form Ai + A2/m (solid line)and i?iln(m) + B2 (dots) with 
selected points of the data in Fig. ^ for the HM (boxes). 
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FIG. 13. Values of Dn for yi up to order 5 in /3 (filled circle) and 11 (filled boxes), and for yi 
up to order 10 in di (empty circles) and with first order corrections in y2 (dots), and up to order 
80 in di (empty stars) and with first order corrections in y2 (solid line). 
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